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The field-cooled (FC) process on an Ising spin-glass model is investigated by a standard Monte 
Cairo (MC) method on one hand, and the equilibrium magnetization of the same system is 
evaluated by the exchang MC method on the other hand. The two types of simulation reveal 
intriguing glassy (nonequilibrium) dynamic properties of the FC magnetization (FCM) of the 
system. Particularly, the FCM decrease is observed when the FC process is halted at a low 
temperater, although its value is smaller than the corresponding equilibrium value. We present 
a comprehensive interpretation of such peculiar phenomena based on the scenario for the FCM 
process recently propsed by Jonsson and one of us (HT). 
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In the spin-glass (SG) study, one of most fundamental 
problems yet unsettled concerns the stability of the equi- 
librium SG phase under a finite magnetic field h. 1 ^ The 
SG mean-field theory predicts its stability up to a certain 
critical magnitude of h, 2 ^ while by the droplet theory 3 ' 4 ' 
it is unstable even in an infinitesimal h. There have ap- 
peared two representative experiments on an Ising spin 
glass in relatively large fields; one 5 -' favors the mean- 
field picture and the other 6 ) does the droplet picture. 
Recently there have appeared numerical analyses on the 
Ising Edwards- Anderson (EA) model which strongly sup- 
port the droplet picture. 7, 8 - ) There still remains, however, 
a most fundamental SG property to be appropriately un- 
derstood, namely, the behavior of the field-cooled mag- 
netization (FCM) under a relatively small field. 

Since the early stage of the SG study, a kink-like be- 
havior in the temperature dependence of FCM has been 
observed, 9 - 1 and it has been considered to be evidence 
for the SG phase transition in a finite field. In detailed 
measurement on a AgFe spin glass, however, Lundgren 
et al 10 ^ found the following peculiar FCM characteristics 
in small h. The FCM exhibits not a cusp but a peak at 
a temperature, denoted here as T*, which is close to T c , 
a phase transition temperature under h = 0. Further- 
more, when the FC process is halted at a temperature, 
denoted by Thit, lower than T* , the magnetization which 
we call the halted FCM (HFCM) initially decreases. By 
later experiments, the HFCM is checked to increase when 
Ti rr < Xhit < T* with T; rr being the so-called irreversibil- 
ity temperature where the zero-field-cooled magnetiza- 
tion (ZFCM) starts to deviate from the FCM. 11 ) The 
changes in FCM involved in these phenomena are small 
in magnitude as compared with a value of the FCM itself, 
but they are considered to be intrinsic properties com- 
monly shared by most of typical spin glasses. Having fur- 
ther confirmed these intriguing FCM properties by both 
experiments on an Ising spin glass Feo.55Mno.4sTi03 and 
by numerical simulations on the Ising EA model, Jonsson 
and Takayama (JT) 12 ' have proposed, based on the SG 
droplet picture, a novel scenario which explains these 



FCM properties as the consequence of the glassy (non- 
equilibrium) nature peculiar to the FCM dynamics in- 
volved. 13 ' 

In the JT analysis, however, important information 
has been missing, i.e., the one on the equilibrium mag- 
netization (EQM) for a given set of temperature and 
field, M eq (T,h), which is the destination of the HFCM 
at the corresponding (T, h). To our knowledge, there has 
been no real experiment, in which the HFCM is observed 
to reach the EQM. This is consistent with the intrinsic 
SG nature, i.e., the marginal stability predicted by both 
mean-field and droplet theories. It implies that relaxation 
times are distributed continuously up to infinity not only 
at the transition point but also in a whole SG phase. By 
computational experiments, on the other hand, we can 
evaluate the EQM of a finite system by using the sophis- 
ticated numerical methods such as the exchange Monte 
Carlo (MC) method 14 ' which artificially accelerate equi- 
libration in SG systems. We can, of course, analyze the 
FCM of the same system by a standard MC method. The 
information on the EQM is quite important particularly 
for analyzing dynamical processes in such a system with 
the marginal stability. In fact, the main result of the 
present work is that the HFCM at T ut < T* does de- 
crease though it is smaller than the corresponding EQM, 
indicating that the phenomenon is far from equilibrium. 
The purpose of the present Letter is to report the FCM 
behavior simulated further in details as well as the EQM 
data, and then to argue comprehensively such intriguing 
FCM phenomena as the FCM peak, the HFCM decrease 
at T h it < T* , and its increase at T* < T h it < T ilT , based 
on the JT scenario. 

The system we examine is the 3-dimensional Ising EA 
model whose nearest neighbor interactions obey a Gaus- 
sian distribution with mean zero and variance J which we 
use as the units of energy and temperature (with &b = 1). 
Its T c has recently been estimated as T c ~ 0.96. 15 ' The 
strength of field h is represented by its associated Zccman 
energy. We adopt a standard heat-bath MC method to 
simulate FC processes, and the time is measured in unit 
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Fig. 1. 



of MC steps per spin (mcs) . A system is cooled by a step 
of AT = 0.01, while it is held by a staying time, t s ty = n 
mcs, at each temperature. We call this a rate n, or rn, 
cooling process. To specify a value of the FCM per spin, 
we take an average over at most last 30 mcs at each T, 
which is much shorter than n. For the FCM analysis, 
systems only with the linear dimension L = 24 are sim- 
ulated as in our previous works. 12 ' 16 ' We take averages 
over several thousands or more samples to get sufficient 
accuracy for discussion of small differences in FCM and 
HFCM of interest. To evaluate the EQM we carry out the 
exchange MC simulation 14 ) with 64 temperature points 
in the range 0.5 <T< 2.75, and take average over 1500 
samples. The periodic boundary condition is imposed for 
systems simulated by both methods. 

Let us begin our discussion with the EQM results 
shown in Fig. 1. The EQM's of the systems with L = 8 
and 24 coincide with each other even at lower tempera- 
tures than T c when h is large (= 0.4), but they don't for 
small h such as h = 0.1. It implies the presence of a cer- 
tain characteristic length scale which is smaller (larger) 
than L = 8 under h = 0.4 (0.1). In the SG droplet theory, 
we immediately think of such a length scale, namely, the 
field crossover length, Lh- It separates the characteris- 
tic behavior of droplet excitations in the equilibrium SG 
state by their size I; it is dominated by the Zeeman en- 
ergy (~ hl d t 2 ) for I > Lh and by the SG stiffness energy 
(~ Yy^ZEA' ) f° r I < Lh. Here, d is the spatial dimen- 
sion, T the stiffness constant of the SG ordering, <7ea 
the EA order parameter, and 9 the stiffness exponent. 
Explicitly, Lh is written as 

-5 



Lh 



IrhT 



(1) 



where <5 = (d/2 — Q)~ l and It is a constant weakly 
depending on T through c?ea and T. Then the L- 
independent EQM at sufficiently low temperatures sim- 
ply implies that Lh < L = 8 for a large field (h = 0.4). 

The simulated L — 24 EQM seen in Fig. 1 is nearly 
constant up to a temperature much higher than T c . This 
observation leads us to extend the presence of Lh to a 
temperature region higher than T c in the following sense. 
As the temperature is decreased from above, the SG 



short-range-order (SRO) grows but is hindered by field 
h around a temperature, denoted as T crs , where the SG 
coherence length, £ C (T) oc (T — T c )~ v , becomes com- 
parable with the extended Lh- These SRO's grow inho- 
mogencously in space due to the randomness in degrees 
of frustration in the SG system, and they are consid- 
ered to behave as randomly frozen clusters at T < T CIS . 
Namely, each such cluster is fluctuating as a whole in 
equilibrium but with a reduced magnetization per spin 
and so the response of the system to h is reduced. The 
consequence is the EQM smoothed-kink behavior seen in 
Fig. 1, from which we roughly specify a value of T crs . The 
parameters T and q^A involved in the extended Lh arc 
considered to be the corresponding quantities associated 
with the SG SRO. Without going into its details, we sim- 
ply assume here that the temperature dependence of the 
extended Lh is significantly weaker than that of £ C (T), 
and that the extended Lh smoothly changes to the orig- 
inal one at low temperatures. It is also noted that the 
above argument is appropriate only for the case L > Lh- 
For systems with L < Lh, we have to take into account 
the finite-size effect introduced by the periodic boundary 
condition. The L = 8 EQM in h = 0.1, which deviates 
from the L — 24 EQM around T crs , is the case. Its details 
will be discussed elsewhere. 

Now let us move to discussion on the FCM behav- 
ior. In Fig. 2 we show FCM's observed in three FC pro- 
cesses with different cooling rates as well as the EQM in 
h = 0.1. Also shown are the ZFCM observed in the rlO 4 
process. One can see that T; rr ~ 1.3 for the rlO 4 process, 
whereas it is out of the figure, T; rr ~ 1.6, for the rlO 2 pro- 
cess (see Fig. 1 in [12]). The important observations here 
are that the FCM's with slow cooling (rlO 4 and r33333) 
exhibit a peak, and that these FCM'S do not exceed the 
EQM even at T ~ T*. The peak is the sharper, the slower 
is the cooling. In contrast to real experiments, 10,12 ' the 
crossing of FCM's below T* is not ascertained in the 
present simulation. The reason can be attributed to the 
field strength, namely h = 0.1 is still large for the cross- 
ing to be observed. The FCM with rapid cooling (rlO 2 ), 
on the other hand, does not exhibit a peak, and it even 
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Fig. 2. The FCM, ZFCM and EQM in h = 0.1 (L = 24). 
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Fig. 3. The HFCM, QUFM and QUPM in h = 0.1 at three T h it's. 
The horizontal line represents the EQM with its error bar indi- 
cated at the right end. 



exceeds the EQM at low temperatures. 

The HFCM results at three Twit's after the rapid (rlO 2 ) 
and slow (rlO 4 ) FC processes are shown in Fig. 3. Also 
shown there are time evolutions of magnetizations after 
an instantaneous quench from the paramagnetic (com- 
pletely random) and ferromagnetic (completely polar- 
ized) states at the same temperature, denoted respec- 
tively as QUPM and QUFM. We sec in Fig. 2 that 
the equilibrium is realized by the slow cooling down to 
T < 1.1 (< T; rr ) within the present numerical accuracy . 
Correspondingly its HFCM at T^it = 1.2 stays constant 
as seen in Fig. 3(a). The FCM in the rapid cooling, on the 



other hand, is significantly smaller than the EQM already 
at T = 1.2 and its HFCM at Thit =1.2 starts to increase 
at time t > n (= 10 2 ). Its equilibration is accomplished 
at t ~ 10 5 , where both QUPM and QUFM also reach to 
the EQM. Interestingly, the QUFM does so from below, 
after it overshoots the EQM. At a low temperature such 
as T = 0.9, the FCM in the slow cooling, having passed 
its peak, is smaller than the EQM (Fig. 2). Still, as seen 
in Fig. 3(c), its HFCM at T Ut = 0.9 definitely starts to 
decrease, which is the opposite direction to the EQM. 
Also the HFCM at T hlt = 0.9 after the rapid cooling 
starts to decrease from its initial value larger than the 
EQM, continues to decrease, and overshoots the EQM. 
We can see also in the figure that the QUFM exhibits a 
minimum at t ~ 10 5 , but it is still significantly smaller 
than the two HFCM's even at t ~ 10 6 , where the QUPM 
is further smaller than the QUFM. These data suggest 
that it will take 10 10 mcs or much more for these mag- 
netizations to reach the EQM at this temperature. At 
Thit = 1-05 shown in Fig. 3(b), a little lower than T* of 
the slow FC process, the corresponding HFCM exhibits 
monotonous decrease. The HFCM after the rapid cool- 
ing, on the other hand, exhibits a subtle upturn after its 
initial decrease. These results strongly suggest that a cer- 
tain qualitative change in nature of the FCM dynamics 
occurs around this temperature also in the rapid cooling, 
though its FCM does not exhibit a peak. 

Now let us interpret the peculiar FCM phenomena so 
far described along the line of the JT scenario. It has 
been well accepted that, at T ~ Tj rr , spin clusters of 
a mean size denoted as £*, go into (release from) their 
thermally blocked state in the cooling (reheating) pro- 
cess of the FCM (ZFCM) measurement. Here 'thermally 
blocked' means that clusters are not fluctuating within 
time scale of t s ty of the FC and ZFC processes adopted. 
Therefore, Tin- depends on both h and i s ty, while T crs for 
the EQM docs only on h. Correspondingly, £* also de- 
pends on h and t sty . Furthermore, we consider it plausible 
that these thermally blocked clusters appear inhomoge- 
neously in space and are isolated with each others. As 
the FC process is continued below Ti rr , clusters whose 
sizes are smaller than £* become thermally blocked. But 
when the FC process is halted by a much longer time 
than t s ty, the clusters already thermally blocked become 
able to fluctuate and the HFCM increases toward the 
EQM. 

According to the JT scenario, the spin clusters so far 
thermally blocked and isolated from each others become 
in touch with at T ~ T*. Then the SG stiffness energy, 
so far effective only within each clusters, now becomes 
effective also between the clusters, and the latter rear- 
range themselves to a configuration which gains the SG 
stiffness energy. This rearrangement is initiated by in- 
teractions with neighboring spins. It therefore starts to 
occur in a length scale much smaller than £*, which in 
turn is smaller than in the FC process in small h and 
with small t sty . Its consequence is, of course, decrease in 
the total magnetization of the system. Even at T < T*, 
however, there may be still smaller clusters which are 
isolated from others and which are going to thermally 
blocked preferentially to the field direction by the tern- 
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perature decrease of AT. In the rapid FC process, the 
latter dominates the former, giving rise to only decrease 
of the FCM increasing rate but not the FCM peak. The 
opposite is the case in the slow FC process, yielding the 
FCM peak. 

As the temperature is further decreased, cluster flips 
governed by the thermally activated process in this tem- 
perature range and below become less and less frequent 
and the FCM becomes nearly constant. But there cer- 
tainly still exist isolated spins (and small clusters) whose 
internal field acting on them, denoted as faint) is van- 
ishingly small. Actually, the distribution of faint's of the 
mean-field model is shown to be an even function of hint 
and is zero at h lnt = in the limit fa —> 0. 18 -* This is natu- 
rally expected to faint 's in the EA model of the present in- 
terest as well. These isolated spins and small clusters are 
expected to be blocked equally to both directions rela- 
tively to the field on average, yielding the constant FCM. 
Similar FCM behavior is observed in a superspin glass 
consisting of magnetic fine particles interacting with each 
others via dipolc-dipolc interactions. 19 -' 

When the FC process is halted at T < T* , the rear- 
rangement of clusters mentioned above proceeds and so 
the HFCM decreases. As seen in Fig. 3(c), this happens 
even at a halt of the rapid FC process after an interval 
shorter than i sty where the thermal blocking of small 
clusters due to the temperature decrease of AT have 
completed. This HFCM decrease is expected to continue 
until the mean size of the SG SRO proper to T = Thit 
and fa = 0, denoted as £r(i), reaches to L] % . At further 
longer time, £r(i) no longer grows. Instead, the system 
consisting of clusters of a mean size Lh relaxes to its equi- 
librium configuration, yielding the HFCM upturn and its 
increase toward the EQM. The HFCM upturn after the 
initial decreases is vaguely seen in the rapid FC process 
as shown in Fig. 3(b). The QUFM upturn after the initial 
rapid decrease seen in Fig. 3 can be similarly interpreted. 

The computational results so far described agree qual- 
itatively with those observed in real experiments on var- 
ious typical spin glasses, though the scales of time and 
field strength between the two types of observation quan- 
titatively differ very much. In our previous work on the 
field-shift aging protocol on Ising spin glasses, 7 ) the com- 
putational result is extended to the laboratory time range 
by making use of the growth law of the above-mentioned 
£r(t) which is also obtained numerically. 16 ) This repro- 
duces even semi-quantitatively the corresponding exper- 
imental result, 5 ) implying that the two measurements 
look at the same physics. We believe that this is also the 
case for the present FCM analysis in a relatively weak 
field. 

Recently non-equilibrium phenomena below T* (~ 
T c ), such as aging, memory, and rejuvenation, have been 
extensively studied. ^ Although most of the results are 
interpreted by regarding that the state reached by a cool- 
ing is in near equilibrium, that is in fact not the case as 
clearly demonstrated by the HFCM decrease at Thit < T* 
in the present work. We consider that quasi-equilibrium 
interpretation is appropriate only for phenomena within 
a time scale of the order of the waiting time adopted, 
such as t s ty in the present analysis. If the measurement 



is continued to much longer time, the non-equilibrium 
nature discussed in the present work is expected to be 
observed. In fact, the post-aging decay phenomenon in 
such a time range has been recently reported. 20 ) We con- 
sider it to be due to the presence and very slow decay of 
clusters at T < T c which have been thermally blocked at 
T ~ Tj rr > T c , or more generally, at T where the char- 
acteristic relaxation time of the SR SRO of a size £ C (T) 
exceeds t s t y of a cooling process adopted. 

To conclude, by computational experiments, or by 
numerically looking at physical properties that a well- 
defined, rather simple theoretical (Ising EA SG) model 
exhibits, we have observed its intriguing FCM phenom- 
ena such as the FCM peak and the HFCM increase (de- 
crease) at a halt at T* < T hlt < T irr (T ut < T*). The 
computational experiment can further evaluate equilib- 
rium properties by artificially accelerating the relaxation 
dynamics. Combined the EQM result thus obtained with 
the FCM behavior, we have presented a comprehensive 
interpretation of most of the properties associated with 
the FC process based on the JT scenario. Although we 
have introduced certain assumptions such as the presence 
of the extended field crossover length up to temperatures 
above T c , we rather consider the latter one of the conclu- 
sions reached by the present computational experiment. 
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